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ABSTRACT 

We use general arguments to show that a continuously powered radiative 
blast wave can behave self similarly if the energy injection and radiation 
mechanisms are self similar. In that case, the power-law indices of the blast wave 
evolution are set by only one of the two constituent physical mechanisms. If the 
luminosity of the energy source drops fast enough, the radiation mechanisms 
set the power-law indices, otherwise, they are set by the behavior of the energy 
source itself. We obtain self similar solutions for the Newtonian and the 
ultra-relativistic limits. Both limits behave self similarly if we assume that the 
central source supplies energy in the form of a hot wind, and that the radiative 
mechanism is the semi-radiative mechanism of Cohen, Piran & Sari (1998). 

We calculate the instantaneous radiative efficiencies for both limits and hnd 
that a relativistic blast wave has a higher efficiency than a Newtonian one. 

The instantaneous radiative efficiency depends strongly on the hydrodynamics 
and cannot be approximated by an estimate of local microscopic radiative 
efficiencies, since a fraction of the injected energy is deposited in shocked matter. 
These solutions can be used to calculate Gamma Ray Bursts afterglows, for 
cases in which the energy is not supplied instantaneously. 


Subject headings: Gamma raysibursts — hydrodynamics — relativity — shock 
waves 


Introduction 


Afterglows from Gamma-Ray Bursts (GRB) have been discovered for 13 GRBs 
since the hrst detection of the afterglow of GRB970228 ( e.g. |in h /and et. al. 1998|) . 
The simplest cosmological hreball afterglow model ([Paczyhski fc Rhoads 19^ , [Katz 
1994| , pVleszaros fc Rees 1997] , jWaxman 1997| , [Meszaros, Rees, fc Wijers 1997] , [Piran 199BD 


seems to be in a good general agreement with the observed behavior ( Wijers, Rees and 


Meszaros 1997[) , or at least with the power-law decay. However, several works have tried to 
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investigate more subtle effects which can change the afterglow characteristics. In particular, 
Rees fc Meszaros (1998)| and [Panaitescu, Meszaros and Rees (1998)| have tried to explain 
the deviations from an ideal power-law, and the variety of light curves using a model in 
which a source emits shells with different Lorentz factors, which results in a gradual energy 
supply to the shell of shocked matter. In their work they have neglected the thickness of the 
shell of shocked matter, and did not calculate its structure. However, as evident from the 
classical treatment of adiabatic blast waves by [Blandford fc McKee (1976]i the evolution of 
a blast wave depends strongly on the detailed structure of the matter inside it. Therefore, it 
remains to investigate the complete hydrodynamics and energy budget of a slow power-law 
decay blast wave. 


The dynamics of blast waves with gradual energy supply is also relevant for the study 
of compact steep-spectrum objects (CSS), Gigahertz peak spectrum objects (GPS), and 
active galactic nuclei (AGN) radio lobes. The generally accepted model for these objects 
(see e.g. pcheuer 1974|, [Begelman fc Cioffi 1989|, [Bicknell, Dopita fc O’Dea 1997|) describes 
a blast wave continuously powered by a jet from an AGN. Similar models have also been 
used for interstellar bubbles ( [Castor, McCray and Weaver 1975|JWeaver et. al. 1977|) , for 
plerionic supernova remnants (see e.g. [Weller 1983|) and for galactic super-shells ([McCray 

[ 19871 ) . 


In this paper we consider blast waves with a gradual energy supply by a hot wind. 
Specihcally we obtain a new self similar solution for a continuously powered radiative 
blast waves. Even if this situation may not be directly applicable to astrophysical objects, 
it is still important. Self similar solutions are simple enough that they may be solved 
analytically (sometimes), and are easy to grasp. Furthermore, it is likely that a generic 
blast wave will have tendency to approach self similarity. 


Adiabatic self similar blast waves with energy injection have been treated previously 
by Blandford fc McKee (197^ in the relativistic regime. In the Newtonian regime Castor, 
McCray and Weaver (1975) and [Ostriker fc McKee (1988]] have investigated both radiative 
and adiabatic blast waves with energy injection, but have not treated the case of fast 
decaying sources, which did not £t into their self similar framework. In this paper we show 
that the inclusion of radiative losses has a qualitative influence on the solutions. Using 
our new solutions we obtain an accurate result for the fully radiative steady injected blast 
waves, which were treated approximately by Blandford fc iVlcKee (197U] . We add new 
solutions for a region of the parameter space which has not been treated yet neither in the 
Newtonian nor in the ultra-relativistic limits. 


Self similarity appears only if the physical situation can be characterized by a minimal 
number of dimensional parameters, such that physical scales can be constructed only from a 
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single combination of those parameters. Therefore, by assuming that a blast wave behaves 
self similarly, we can deduce the functional form of different physical processes. In fact, 
adiabatic blast waves already exhibit self similar behavior without additional parameters. 
This requires that energy injection and radiation processes in self similar blast waves will 
not have intrinsic scales. Specihcally, the injection rate must be a power-law. 


The radiation process should also be self-similar. In order to satisfy that demand, we 
assume that the radiation process results with a semi-radiative scenario where the cooling 
mechanisms are fast comparable to hydrodynamic time scales (“fast cooling”), but only 
a fraction of the energy produced by the shock is radiated away. This would take place, 
naturally, in any collisionless shock acceleration (see |(Johen, Fir an fc Sari 199^ , hereafter 
CPS98). For example, in a GRB afterglow the cooling time scales are shorter then the 
hydrodynamic time scale, (see e.g. |Waxman 1997|, [IVleszaros, Rees, fc Wijers 1997|, |Sari 


Piran, fc Narayan 1998 ), and a fraction of the energy stays with shocked protons which do 
not cool. 


We describe our model in Sec. It is composed of a central source emitting a wind 
which interacts with a surrounding medium and creates a blast wave. We proceed in Sec. 
^ by calculating the blast wave energy and the radiated luminosity assuming self-similarity 
alone. In Sec. ^ and Sec. we split the discussion to the Newtonian and ultra-relativistic 
limits, assuming the self-similar radiation mechanism of CPS98. For each case we obtain 
the hydrodynamic solution and the radiative efficiencies (analytic solution for the ultra 
relativistic limit, and numerical for the Newtonian limit). Finally we summarize our results 
in Sec. ||. 


2. The Model 

We consider a spherical semi-radiative blast wave, that appears when energy is released 
continuously into an ambient medium. This results in a strong shock wave that expands 
supersonically. We consider the regime where the influence of the injected and initial mass 
is negligible, and that the pressure of the surrounding medium is small compared to the 
energy density of the flow. These assumptions are necessary in order to obtain a self similar 
solution (see e.g. 

We assume that the source supplies energy in the form of a hot wind with a negligible 
mass. This wind pushes away the surrounding medium, and creates a cavity. This simple 
model leads to different qualitative behaviors, depending on the expansion velocity of the 
blast wave. If this velocity is much lower than the speed of light, (Newtonian limit) the 


Landau fc Lifshitz, §99). 
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sound velocity of the ejected wind can be much higher than the expansion velocity. In this 
case the pressure inside the cavity settles fast and becomes isobaric. This isobaric bubble 
pushes away the ambient matter and creates a shock wave. (This model, where the injected 
wind is treated as an isobaric interior has been widely used in the study of radio lobes in 
active galactic nuclei (^cheuer 1974|) ). 

In the ultra relativistic case the hot wind has different dynamics. The sound velocity of 
the wind is of the order of the speed of light, but it is always comparable to the expansion 
velocity. In this case we assume that the wind is so sparse that its’ particles inside the 
cavity do not interact, and simply move in straight lines, close to the speed of light, until 
they reach the edge of the cavity. Near this edge they undergo a strong shock, presumably 
a collisionless one, which compresses them. Consequently, this compressed fluid produces a 
shock which advances into the ISM. 

Both in the Newtonian and in the ultra relativistic limits the blast wave is composed 
of four regions, (a) The injected wind: contains wind that did not interact, yet, with the 
blast wave. In the Newtonian limit the reverse shock reaches the center, and this region 
disappears, (b) Shocked wind: contains wind that was compressed due to the interaction 
with shocked ISM. It is separated from region (a) by a reverse shock, (c) Shocked ISM: 
contains ISM that has been compressed by the forward shock at the blast wave front. A 
contact discontinuity separates this region from the shocked wind region, (d) ISM at rest. 
A schematic view of the four regions appears in Fig. |I|. 


3. Self similarity in energy injection cases 

We look for self similar solutions for semi-radiative blast waves with central energy 
supply. If a blast wave contains only one of these two mechanisms ( energy injection and 
radiation ), each of them leads to a different temporal behavior. As self similarity does not 
allow the existence of several time scales in the solution, it is apriori not known whether 
such solutions exist. 

Self similarity requires a dimensionless energy injection mechanism. We therefore 
deduce that the energy is supplied to any interval along the self similar solution at a rate 
L (X t^. The energy transfer rates along the self similar hydrodynamic prohle ( especially 
Lrs and Led ) scale the same, and differ only by a constant factor, which depends on 
the exact solution. However, this scaling can differ from L^, the energy supplied by the 
central source, as a fraction of the injected energy resides in the interior regions. Region 
(a) contains a freely expanding wind. If it exists, it does not depend on the blast wave 
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parameters. For self similarity the energy transfer rate from region (a) to region (b-c) must 
satisfy the self similarity scaling laws. However, it is possible that the energy in region (a) 
will evolve differently. 

Assuming that a self similar solution exists, we can write the energy loss rate as 

E 

Lrad '' ^ ^ b) (1) 

where E is the energy stored in a hxed interval of the self similar prohle. In the Newtonian 
limit we choose E to denote the energy stored in shocked ISM alone (and the corresponding 
injection rate is L^d) and in the Ultra-relativistic case E denotes the energy stored in the 
entire shocked matter (and the corresponding injection rate is L^s)- The value of the 
constant k depends on this dehnition, but it is always dimensionless. This proportionality 
is evident from self similarity, as there exists only one way to construct the dimensionality 
of Lrad (energy/time) using the characteristic parameters. Note that this argument holds 
for any self similar radiation mechanism. 


The energy conservation equation includes losses as well as injection from the blast 
wave center. 


(W 

dt 


lLy 


E 

K —. 

t 


( 2 ) 


For a ^ —1 — K, Eq. ^ has an analytic solution 


E 


^ -i^rt+At- 


K + a + 1 tE 


( 3 ) 


where A is set by the initial conditions. 


We proceed by investigating the asymptotic behavior of this equation for different 
values of a and k. In cases where the equations behaves self similarly we hnd the power-law 
index of the energy temporal behavior E oc t^, and obtain the relations between the various 
power-law indices. A summary of the different limits appears in table 1. 

For L = 0, we reproduce the instantaneous energy injection case and Eq. trivially 
becomes a power-law with A = —k. Even in this simple case, an extrapolation toward 
f > 0 results in an inhnite energy. From Eq. |l| it is obvious that the radiated energy also 
diverges in this limit. These two inhnities cancel each other, and at hnite times the energy 
is always hnite. A specihc blast wave with an arbitrary initial energy evolves to this self 
similar behavior at late times. Clearly the limit f 0 has no physical relevance and only 
the late time behavior is interesting. It is therefore better to treat the solution from t = oo 
backward in time, as this is the region where we expect the blast wave to evolve according 
to the self similar solution. At any hnite time all the energies are hnite, and we encounter 



no infinities using this view. We discuss which is the energy that will be radiated 

from time t to inhnity. It is evident that E^^{t) = E{t), as the energy stored in the shell 
decreases to zero when t ^ oo. 


Sources which scales with a < — 1 were not treated so far. In this scenario the injected 
and radiated energies from f = 0 are inhnite, similarly to the L = 0 case. This means that 
the solution is not physical near t = 0. As in the L = 0 case, we overcome the inhnities 
by discussing E'^-{t). The solutions are divided into three sub-classes, depending on the 
relation between a and k\ 

• \i a > —1 — K the injected energy decreases slower than the rate of an un-injected 
blast wave, and the power-laws are set by the injection law [|. The energies injected 
and radiated from t = 0 are inhnite. However, the energy that will be radiated until 
inhnity is always larger than the stored energy, and the energy which will be injected 
is larger (smaller) than the stored energy if a is larger (smaller) than —1 — k/2. Note 
that even if the stored energy is larger than the energy injected until inhnity, the 
evolution is set by the injected energy power-law. In both cases A = cr -|- 1. 

• a = —1 — K then Eq. is not valid. The solution of Eq. ^ for this case is 



( 4 ) 


This equation has no power-law asymptotics, and there is no self similar solution in 
this case. 

• In cases where a < —1 — k, the late time behavior is identical to the instantaneous 
injection case, and the asymptotic of Eq. ^ is a power-law with A = —k. Note that in 
all the sub-cases with a < — 1 the stored energy drops to zero with time. 

Sources with a > —1 were treated by [Ostriker fc McKee (1988)i In this case 
the energy injected from f = 0, is hnite for every t, but it diverges as t —> oo. This 
divergence allow us to neglect any initial energy in the blast wave. Eq. ^ again has a 
power-law asymptotic, now with A = cr -|- 1. The energy stored in the shocked matter is a 
constant fraction of the energy injected so far, and the same holds for the radiated energy. 
The energy injected is always larger then the energy stored in the shocked matter, and the 
radiated energy from t = 0 is larger (smaller) than the stored energy depending if k is larger 
(smaller) than a -|- 1. 


^ Even though the power-law index is set by the injection mechanism and not by the radiative mechanism, 
the radiated energy is not negligible. 
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If the energy is supplied with cr = — 1, the integral of injected energy is logarithmic, 
and it diverges in both limits of t. However, the asymptotic of Eq. ^ is again a power-law, 
now with A = 0, i.e. the shell holds a constant energy which is set by initial conditions. 

The study of the energy conservation equations revealed the temporal behavior of the 
energy if the solution is self-similar. However, k cannot be determined by self-similarity 
arguments alone. It is generally determined only by solving the complete hydrodynamic 
equations. We split the discussion of this solution to the Newtonian and the ultra-relativistic 
limits. 


4. The Newtonian solntion 

To solve the hydrodynamic equations we use our model (Sec. |^), where the energy from 
the central source is supplied by an internal fluid with a high temperature and a low density, 
in the limit of inhnite sound velocity. The large sound speed of the injected matter prevents 
the existence of an interior shock, and the flow contains only three regions (see Fig. ^): (b) 
an isobaric interior composed of shocked wind; (c) shocked interstellar medium; and (d) 
ambient ISM. The hot internal fluid pushes the shocked matter at the contact discontinuity 
and accelerates it. The shocked matter accretes mass from the ambient medium, and heats 
it at the shock front. The heated matter emits a fraction of its internal energy near the 
shock front. 

If the evolution is governed by the injection mechanism, the blast wave is characterized 
by its instantaneous radiative efficiency, which is the ratio between emitted energy (at the 
shock front) and the energy released at the center, at the same time. Due to the structure 
of the blast wave (the three regions), it is helpful to divide the energy transfer into two 
stages: (i) energy transfer from the wind to the shocked matter; and (ii) energy transfer 
from the shocked matter to radiation. 

Both stages show instantaneous efficiencies which are larger than unity if the source 
luminosity drops fast enough. This is caused by the delay between the time energy is 
supplied to the shocked matter and the time this energy is radiated. For example, if the 
source shuts down spontaneously it emits no energy, but energy is still radiated from the 
hot shocked matter. This results in an inhnite instantaneous efficiency. 

During the expansion, the hot internal huid (region (b)) performs work on the shocked 
matter (region (c)). Led = ^T^fldPcdVcd, where rcd,Pcd, and Ved are respectively the radius, 
pressure and velocity of the contact discontinuity. Meanwhile, the internal huid absorbs 


energy from the internal source which emits energy at a rate 


Lin = 47r/3[ 


d{rldPcd/{l - 1 )) 

dt 


drld. 




( 5 ) 


Prior to solving the self similar hydrodynamics, we use the self similar relations to calculate 
the efficiency Led/Lin. A Newtonian blast wave must evolve according to the self similar 
variable (see e.g. [Landau fc Lifshitz| , §106): 




Pi 


1/5 




/ Rsfi, 


( 6 ) 


where pi is the density of the surrounding medium, and Rgh is the radius of the shell. 
The self similar quantities, which turn the hydrodynamic equations into coupled ordinary 
differential equations (ODEs) (see e.g. [Landau fc Lifshit^ , §106) are 


p{r, t) = 

(^) pMO, 


\ 7 -lJ 


u{r,t) = 

2 r 

- 1 1 ^•9^ p 

7 -Fl Rsh{t) 

( 7 ) 

p{r,t) = 

7 -Fl Rsh{t) 



Keeping in mind that the ratio Ved/Rsh is constant (again due to self similarity), we 
substitute p(r, t) into Eq. || and use Eq. ^ to obtain 




Led/L 


in 


3(7-1)(A + 2) 
2(A-3) + 37(A + 2)' 


( 8 ) 


For A > 0, (which corresponds to energy injection cases with a > —1), rjin < 1, as expected. 
However, for A < 0, > 1, as the radiated energy is mainly supported by energy that is 

already in the shocked matter. For A = — 6(7 — l)/(2 + 37 ), pin is infinite. This corresponds 
to the [Ostriker fc McKee (1988)[ solution of adiabatic interior with no internal energy supply 
which pushes a fully radiative shell. 


This efficiency reveals the ratio between the injected energy Lin and L^d, the rate in 
which energy is supplied to the shocked ISM. To obtain the instantaneous global efficiency 
we use Eq. and find that rjrad = Lrad/Lcd = h/{k + cr + 1). 

To obtain k we need the full hydrodynamic solution. We use the energy loss rate of the 
semi-radiative model (CPS98) 


Lrad — ‘^'^PlUshd^sh^i 


( 9 ) 
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substitute the self similar relations of Eq. |^, and obtain 

K = -27r(^-|^)^eoe- 
The instantaneous efficiency is than 

_ Lrad _ 3(7 — 1)(A + 2)/[2(A — 3) + 37(A + 2)] 

~T~~ 1 +125(1+ a)/[(3 + a)327reo^e] 

where ^0 is the self-similar position of the shock. 

To obtain ,^0 and the instantaneous efficiency, we solve the self-similar equations. The 
solution is obtained by integrating numerically the self similar hydrodynamic ODEs from 
the shock down to the contact discontinuity (where the fluid does not move on the self 
similar prohle). The boundary conditions for this integration are the semi-radiative shock 
conditions (CPS98) 

^ '^(^ 0 ) = 1 + — S, p{^o) = 1-1- — S, (12) 

where 6 is related to the radiative efficiency by 

«W = tA^|2+(7-1)'5]. (13) 

1 + 7 

We hnd ^0 using the normalization condition which equates the energy dehned in the 
self-similar parameter (Eq. P) to the energy stored in the shocked matter 

E(t) = (14) 

If the evolution is governed by the radiation mechanism, then A is not known prior to 
solving the hydrodynamics. In that case in addition to Eq. |T^ we use E = —Lrad, and 
obtain both A and ^0 simultaneously. 

We present the hydrodynamic solutions for several cases in Fig. and in Table 
In all cases where the evolution is set by the energy injection mechanism, a contact 
discontinuity appears in the solution. If the injected energy decreases fast enough with 
time, the shocked ISM prohle is similar to the one found in instantaneously injected blast 
waves. The temperature rises monotonously toward the center, and it diverges at the 
contact discontinuity. The velocity is monotonous, and the huid near the shock front 
has the maximal velocity. At the contact discontinuity the huid is at rest relative to the 
prohle. The density is also monotonous, and it reaches the contact discontinuity with a zero 
gradient. This results in a smooth transition between the shocked wind (which composes 


( 10 ) 

( 11 ) 
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the internal isobaric bubble), and the shocked ISM. Effectively, there is no discontinuity at 
the contact discontinuity. 

The boundary between the shocked ISM and the internal hot isobaric bubble depends 
on the energy injection rate. As the injected energy becomes more dominant, i.e., it 
is injected with higher a, the matter tends to concentrate in a narrower shell near the 
front shock. If a < (8 — 97)/(4 + 87 ), the density gradient still vanishes at the contact 
discontinuity, which contains no jumps. For a larger than this value, the density gradient is 
infinite, and the shocked matter region ends abruptly at the contact discontinuity. If cr > 2, 
the density it self also diverges to infinity. Despite these infinities the total mass and energy 
stored in the self similar solution are finite, and the solution is a proper one. Increasing 
a also changes the velocity profile considerably, as for large enough a the maximal fluid 
velocity is not at the shock front, but at the contact discontinuity. 

Inclusion of a radiative mechanism decreases the shell width, with the limit of zero 
width if e —>■ 1. The qualitative effect of radiation is mainly observed in cases where a < — 1. 
There, the solution is identical to the evolution of a blast wave with instantaneous energy 
release, unless the radiative efficiency is high enough. If k > 1 + a, the evolution is set by 
the injection mechanism, and the blast wave has the same structure as ones with a > — 1 . 

The instantaneous efficiency as a function of a is depicted in Fig. The efficiency 
is larger than unity if a < —1 as the main source for radiation is the hot shocked matter 
which was heated at earlier times, and not the energy which is currently being supplied. If 
a > —1, the main contribution to the shocked matter energy is the recently released wind. 
The radiated energy is therefore directly connected to the wind released at that time, and 
the efficiency is lower than unity. It appears that the hydrodynamics are an important 
sink for injected energy. For example, a steady source with a semi radiative mechanism of 
e = 0.5, results in an overall efficiency of 0.34 if 7 = 5/3, and 0.19 is 7 = 4/3. 

The limit of e = 1 cannot be calculated numerically because of the divergence in the 
semi-radiative boundary conditions. However, a complete analytic solution exists for this 
case (see appendix), assuming that the whole ISM is concentrated in a narrow shell. The 
limit of the numeric solution with e —> 1 reaches this solution. 


5. The Ultra-Relativistic solution 

If the blast wave expansion velocity is ultra-relativistic, the sound speed of the wind 
cannot be considered infinite, an internal shock exists, and the blast wave exhibits the 
four regions discussed in Sec. 0. The inner region of the blast wave contains wind which 


11 


flows unaltered inside the cavity until it reaches the blast wave edge. This wind has no 
dependence upon the blast wave parameters, and it can therefore scale differently in time. 


In the derivation of the hydrodynamic solution we follow the footsteps of [Blandfoi^ 

&; McKee (1976)| . We assume that the inner source varies with time as Lin = L^tl, and 
look for a self similar solution for the shocked matter (wind and interstellar medium) only. 
According to our model, the shocked matter flows adiabatically in the entire self similar 
region of the blast wave. The boundary conditions for the solution are found using the 
internal wind, as the energy supplied to the shell by the wind is equal to the sum of the 
radiated energy and the energy stored in the shell. Boundary conditions for the forward 
shock are found using the semi-radiative jump conditions (CPS98). 


Self similar solutions for relativistic blast waves were found ( [Blandford fc McKee 
( 1976)1) to scale according to the self-similar parameter 


X = [1 -I- 2(m -I- l)r^(t)](l — r/t), m > —1 


(15) 


where T is the Lorenz factor of the shock, which scales as T^ cx t~'^. A blast wave which 
evolves with m < — 1 has no time-like relation between the central source and the blast 
wave edge, and thus cannot be reached by injecting energy from a central source. However, 
in other cases, especially with a steep ambient density gradient, a faster acceleration is 
possible. This type of solutions will be discussed in Cohen, Piran & Sari (1999). 

We start by identifying the different regions in the self similar solution. The reverse 
shock, which separates the self similar portion of the blast wave from the inner wind must 
evolve self similarly, and occur at a hxed self similar location y = Xrs- The velocity of 
the reverse shock is therefore v = dr{xrs,t)/dt, and its Lorentz factor is P^s = P/yWI- 
Moreover, due to the high velocity of the wind, the reverse shock is strong and the shocked 
fluid must leave the shock with a velocity of c/3 (|Landau fc Lifshft^ , §135), which translates 
into a Lorenz factor of 7 = v^Prs in the unshocked fluid frame. 

We can therefore use the self similar substitutions for the pressure (p), Lorenz factor 
( 7 ), and density (p) (BM76), 

P{r,t) = ‘^piT{tff{x), 

= ^T{tyg{x), (16) 

p(r,t) 7 (r,t) = 2piT{tfh{x) 


to obtain a condition for the reverse shock, g{Xrs)Xrs = 4. Similarly, for the contact 
discontinuity, where the fluid is at rest relative to the self similar prohle, we obtain 
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g{Xrs)Xrs = 2. These relations enable us to identify the composition of a certain self similar 
interval (shocked ISM or shocked wind) simply by the term g{x)X- 

The wind, which advances with a Lorenz factor which is much larger than T^, overtakes 
the shocked shell with a velocity difference of l/2r^. The rate in which energy is supplied 
to the shell is therefore dE/dt = Lin{tf,)/2Tl, where accounts for the delay between the 
emission of the energy, and the impact at the reverse shock. 

We start by obtaining an equation of motion using the energy conservation equation. 
The energy stored in shocked matter is 

pR{t) 

E= —r^e(47^(r, t) — l)(ir = 87rpiasr^f^/3(m + 1), (17) 

Jo 3 

where we dehne the dimensionless energy integral 



Using this energy and the semi radiative model we write the energy conservation equation 
as 

dE/dt = L{te)/2T‘l — 87ipiRlfJP‘^e/3, (18) 

and obtain 

^-^0 Y/(q+2)^iq-2)/(q+2) ^ 

M67rp ’ 

where 

K = [ 23 '?(a,(g + 1) + e)/{q + 2)'']-V(<?+2)^(.+i)/(.+2)_ 

To hnd as we solve the hydrodynamic self similar ODEs (The analytic solution is given 
in appendix B.). The boundary conditions at the shock front {x = 1), are the semi-radiative 
boundary conditions (CPS98) 

/(1) = 1, J(l) = l + e, = (20) 

The hydrodynamic prohle begins at the shock with these boundary conditions, and 
9{x)x = 1 + e- Near the shock the fluid flows slower than the self similar prohle. Toward 
the blast wave center the huid decelerates relative to that prohle, where at g{x)x = 2 the 
huid is at rest. The contact discontinuity resides in this location. Closer to the center, the 
huid accelerates (the huid here moves faster then the self similar prohle), until it reaches 
the reverse shock where g{Xrs)Xrs = 4. 

We present the hydrodynamic solutions for several cases in Fig. and in table |^. 
The solutions show the same qualitative behavior as the Newtonian blast waves. If the 
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evolution is set by the injection mechanism, the contact discontinuity and the reverse shock 
exist. However, if the injected energy decreases fast enough with time, the shocked ISM 
profile is similar to the one found in blast waves without a continuous power supply. The 
density, pressure and velocity are monotonous, and reach the maximum at the shock front. 
The density vanishes at the contact discontinuity, with a gradient which diverges to —oo. 
The same structure appears also for q < —1, as long as the evolution is set by the injection 
mechanism. If q > 2, the density itself diverges to inhnity at the contact discontinuity. High 
injection rate also results in a velocity profile with a maximum at the contact discontinuity. 


An immediate consequence of the solution applies for the steady injection fully radiative 
limit. [Blandford fc McKee (197^ (Eq. 82 there) have treated this case assuming that 
the entire radiation is emitted from some average location within this layer. We take into 
account the momentum losses due to radiation through the entire cooling layer, and obtain 


r* = 0.79Li„/87rpifi^ 


( 21 ) 


Similarly to the Newtonian case, we would like to characterize a blast wave by its 
radiative efficiency. Using Eq. o and Eq. ^ we calculate the ratio between the radiated 
energy rate and the energy supply rate to the shocked matter (which includes shocked ISM 
and shocked internal fluid), and obtain 


Lrad 1 

Lrs l + as{l + q)/e 

Note that prior to the substitution of as from the full hydrodynamic solution, we can use a 
zero order approximation and assume that as is constant. We can then immediately deduce 
that similarly to the Newtonian case, the efficiency rises with increasing e and decreasing q. 

The instantaneous efficiency is intrinsically ill-defined in the ultra relativistic case. 
The luminosity scales with time as t(2+3g)/(2+q) Eq. |T^), which is generally different 

from the source behavior Therefore, the overall radiation efficiency is not constant in 
time and cannot be used to characterize the blast wave. However, due to the relativistic 
motion of the shell, the radiation it emits is measured by the observer at a different time, 
tabs = t/2/{m + l)/r^. Still, the central source does not move. If it had been measured 
directly by an observer, it would evolve with f, and not tobs- 

We dehne the luminosity in a more observable fashion, by comparing the luminosity 
measured by an observer to the luminosity that would be measured if the blast wave shell 
had not existed, and the source energy could escape freely to the observer. Transforming 
the radiated luminosity to the observer frame using Eq. and dividing by the source 


( 22 ) 
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energy injection rate we obtain: 

_ Lrad _ /"OQ' 

Vtot — — z ^ yj 5 (23^ 

Lin 1 + 0:^(1 + q )/e 

which is constant in time. Efficiencies for different cases appear in Fig. Note that the 
instantaneous efficiencies are higher than those of Eq. ^ which means that energy is 
supplied to the shocked matter faster than it is emitted. This apparent discrepancy is due 
to time contraction, which causes energy emitted over a long duration at the source to 
concentrate into a short observed time. 


6. Discussion and Conclusions 


Radiative blast waves with energy injection exhibit self similar relations both in the 
Newtonian and in the ultra-relativistic limit. We have found that if a self similar solution 
exists, the power-law index can be one of two: (i) The power-law index of a radiative blast 
wave with the same radiative mechanism but without energy injection; (ii) The power law 
index of the central source, disregarding radiative mechanism. The chosen power law is the 
one which corresponds to an evolution with a faster energy increase, or a slower energy 
decrease. Sources which emit energy with a steep decreasing power law (a < —1), behave 
according to the same rule. This corrects a common belief which states that such sources 
can be considered as instantaneous energy release sources. We have shown that this view is 
correct for adiabatic blast waves, but that it does not hold for radiative blast waves. 

The radiative efficiencies in the ultra relativistic case are higher than those of 
Newtonian blast waves. Therefore, the possibility that GRB afterglows are powered by 
continuously emitting sources ([Katz fc Piran 1997| ) cannot by ruled out on efficiency basis. 
Moreover, the instantaneous efficiency of a blast wave with a tuned power law behavior can 
signihcantly exceed the efficiency of the microscopic radiative mechanism. 


A. Analytic solution for the Newtonian fully radiative limit. 


In the fully radiative case the shocked matter is concentrated on a shell of zero 
thickness, where the interior contains hot isobaric matter. In this case the solution can be 
calculated analytically. 

The interior pressure is: 


p _ (7 - 1 ) 


Vit) 





(Al) 
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where V{t) = AnR^f^/S is the blast waves volume, and E^it/to)^ is the energy of the whole 
blast wave. Using the equation of motion 


d{piVUsh) 

dt 




and Eq. we obtain the radiative efficiency 

_ Lrad _ 9(7 — 1)(2 + A)^ 

Lrad + E 2(3 + 4A) [2(A - 3) + 37(A + 2)] 


(A2) 


(A3) 


B. Analytic solutions for a semi-radiative ultra-relativistic blast wave. 


The self similar ODEs for the adiabatic interior of a blast wave are ( Blandtord fc 
McKee (1976)|) : 


Giy) 

Hy) 

Tdiv) 


1 ding 

9 dx 
1 din/ 

9 dx 
1 din h 
9 dx 


(7m + 3/c — 4) — (m + 2)gx 

(m + l)(4-8^X + ^V) 

8(m — 1) + Ak — {m + k — ^)gx 

{m + l){A-%gx +g‘^X^) 

2(9m + 5fc — 8) — 2(5m + Ak — R)gx + {ra + k — 2)g‘^x^ 
(m + 1)(4 - 8gx + - gx) 


(Bl) 


where y = g{x)x- These equations are applicable for external density gradient of p oc r 
Following piandford fc McKee (197^ we use the relation 


We obtain 


g(v)^ = _^_. 

^ dy vS(v) + 1 


din f 

r(v) 

dy 

yG{y) + 1’ 

dlnh 

'Eiy) 

dy 

yQ{y) + 1’ 

dlnx 

1 

dy 

y{yQ{y) + 1) 


(B2) 

(B3) 

(B4) 

(B5) 


Using Eq. [B^ , which does not appear in IBlandford fc McKee (1976)1 , we obtain x{gx) 
and a full solution, parameterize by gx- As gx has physical meaning (Sec. this 
parameterization is even easier to use than the solution as a function of the self similar 
parameter. 
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We integrate Eq. from 1 / = 1 + e, which corresponds to the the blast wave edge, 

toward the contact discontinuity (?/ = 2), and the reverse shock {y = 4). We obtain 


X2{y) = ej^d\nf/dy)dy ^ 

X:i{y) = = 


„.^(16+28m+m^—3fcm)/2/Q n—ml2—l 

Ul P ? 

^(64—24m—m^+fc(2m—32+3fc))/2/a ^(—4—3m+m^+fc(m+l))/2/(m+l) 
^—hilh 2 lo.^y 2^2(fc—m)//i2^(24—13m+m^+fc(—19+4m+3/c))//i2 




(B6) 


using the following dehnitions: 


hi 


a = (160 + 40m + m^ — 6fc(12 + m) + 9/c^)^'^^, 

/3 = —4 — 4m + y{12 + m — 3k) + y^, 

_ a + 3k — m — 2y — 12 
a — + m + 2?/+ 12’ 

2(480 — 368m + 13m^ + m^ + fc(—464 + 136m + m^ — 3/c(—41 + 3m + 3A;))), 

/i 2 = 2(—12 + 3k + m). 


(B7) 

(B8) 

(B9) 

(BIO) 

(Bll) 


Finally, we obtain the hydrodynamic prohles 


X{y) = Xi{y)/Xi{l + e) 
f{y) = X 2 {y)/X 2 {l + e) 

h{y) = Xs{y)/X,{l + e)^^ (B12) 

9{y) = y/x{y)- 


We have used these relation to calculate the prohles in Fig. As discussed in 

Sec. the parameter y distinguishes between the different regions within the blast wave. 
For example, the dimensionless energy (Eq. 0 ) held within shocked matter (between 
1 + e < y < 4) is simply 

as= f f{x)9{x)dx = [ 1 

Ji Ji+el + yG[y) 

without the explicit usage of the self similar location of the reverse shock. 
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Fig. 1.— A schematic drawing of the regions in an energy injected blast wave. The regions 
are: (a) injected hot wind; (b) shocked wind; (c) shocked ISM; and (d) ambient ISM. The 
boundaries between the regions are a reverse shock at Rrs, a contact discontinuity at Red, and 
a forward shock at Rsh- The central source ejects massless hot wind with a luminosity of Lj„. 
Energy is transfered at a rate Lrs to the shocked matter and at L^d to the shocked ISM. The 
blast wave emits radiation with a luminosity Lj-ad from the shock front. In the Newtonian 
solution E denotes the energy stored in shocked ISM, where in the ultra relativistic limit E 
includes the energy stored in shocked wind as well. Region (a) is negligible in the Newtonian 
regime due to the inhnite sound velocity of the wind and Lrs = Lin. 






Fig. 2.— The pressure p behind the shock front (.^ = ^o) normalized to the pressure near 
the shock in an adiabatic blast wave, for four cases: (i) energy injection L oc adiabatic 
(solid); (ii) energy injection L oc adiabatic (dashed); (iii) energy injection L oc 

semi-radiative e = 0.4 (dashed-dotted); (iv) energy injection L oc semi-radiative e = 0.4 
(dashed double-dotted). The blast wave center (.^ = 0) appears on the right side of the graph, 
where as the front shock (.^ = ,^o) is approached on the left. The isobaric region contains the 
shocked wind, where the rest of the blast wave contains shocked ISM. 
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Fig. 3.— The normalized density p behind the shock front = ^q) for fhe cases shown in 
Fig. ||. Cases with density that drops to zero correspond to a < 3. At higher injection rates 
the density diverges to inhnity at the contact discontinuity. The density at the shock front 
(■C = 'Co ) is set by the semi-radiative conditions, and it does not depend on the injection 
mechanism. 
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Fig. 4.— The normalized velocity v behind the shock front (.^ = .^o) for the cases shown in 
Fig. All the graphs coincide on the curve v = Ugh which is relevant both for the 
interior hot bubble and for the contact discontinuity. 








Fig. 5.— The instantaneous radiative efficiency vs. energy injection rate of a Ne' 
blast wave, for various semi-radiative efficiencies: (i) e = 0.1 (dashed double-dott* 
e = 0.5 ( dashed-dotted); (hi) e = 0.9 ( dashed); (iv) e = 1 (solid). 
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Fig. 6.— The normalized pressure / behind the shock front {x = 1) for four cases: (i) 
energy injection L oc semi-radiative e = 0.5 (solid); (ii) energy injection L oc 
adiabatic (dashed); (iii) energy injection L oc semi-radiative e = 0.5 (dashed-dotted); 

(vi) energy injection L oc adiabatic (dashed double-dotted). 
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Fig. 7.— The normalized Lorenz factor g behind the shock front, for the cases shown in 
Fig. 1^. The reverse shock appears at g{x)x = 4. 
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Fig. 8.— The normalized density h behind the shock front, for the cases shown in Fig. ^ 
Cases with density that drops to zero correspond to g < 2. At higher injection rates the 
density diverges to inhnity at the contact discontinuity. 
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Table 1. Characteristic relations in continuously powered blast waves. 


condition 

^stored 

^ rad! ^stored 

^ rad/ ^stored 

Einj/Estored 

E'^j/E stored 

a > —1 

K+a+l\tL^ 

K 

1 + cr 

00 


00 

a = —1 

L f 
K+cr+1 

00 

00 

00 

00 


L ( i Yt 
K+a+YtL^ 

00 

K 

|l+^l ^ ^ 

00 

K 1 

k+i| 

1 

1 






— 1 — K > a 


00 

1 

00 

0 

L = 0 


00 

1 

0 

0 


Table 2. Parameters of similarity solutions for radiative continuously powered Newtonian 

blast waves with 7 = 5/3. 


A 

a 

e 

1 - U/^o 

^in 

^tot 

4 

3 

0 

0.08 



1 

0 

1/2 

0.05 

6/11 

0.34 

1/2 

-1/2 

0 

0.20 



1/2 

-1/2 

2/5 

0.10 

2/3 

0.23 

-1/4 

-5/4 

2/5 

0.62 

14/9 

8.35 

-4/7 

-11/7 

1 

0 

1 

CX) 










Table 3. Parameters of similarity solutions for radiative continuously powered 

ultra-relativistic blast waves. 


m 

q 

e 

Xcd 

Xs 


CXg 

K 

'f]rad 

^tot 

-1 

oo 

0 

1 

1 

0 

0 

oo 

0 

0 

2 

-2/3 

0 

3.37 

6.28 

0.89 

1.24 

7.50 

0 

0 

2 

-2/3 

1 

1 

1.86 

0 

1.20 

2.22 

0.71 

0.88 

3 

-1 

1 

1 

2.38 

0 

1.41 

8 

1 

1 

4 

-6/5 

1/2 

2.59 

8.1 

0.93 

1.42 

127.46 

2.32 

1.53 

7 

-3/2 

1 

1 

9.04 

0 

1.68 

17417.80 

6.20 

2.06 






